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The  dependence  of  self-pumped  photorefractive  two-beam  coupling  in  Fe:KNb03  with  crystal  angle  has 
been  measured  at  532  nm  for  the  a-c  and  b-c  crystal  planes  using  a  combination  of  oil  immersion  and  relay 
imaging.  A  significant  deviation  between  existing  theory  and  experiment  is  found  for  large  crystal  angles.  We 
propose  this  is  due  to  a  highly  anisotropic  effective  charge  trap  density  in  Fe :  KNb03  of  almost  two  orders  of 
magnitude  variation  with  crystal  angle  in  the  a-c  plane. 

DOI:  10.11 03/Phy sRevB . 73 . 1 74 1 02  PACS  number(s):  61.72.Ss,  42.65.Hw,  42.70.Nq,  71.55.— i 


INTRODUCTION 

Nonreciprocal  light  transmission  through  crystals  of  iron 
doped  potassium  niobate  (Fe:KNb03)  occurs  by  photore¬ 
fractive  beam  coupling  between  two  antiparallel  beams 
formed  from  the  incident  beam  and  a  Fresnel  reflection  from 
the  rear  of  the  crystal.  Interference  between  these  two 
counter-propagating  beams  writes  a  reflection  grating  which 
amplifies  one  beam  and  attenuates  the  other.  The  direction  of 
amplification  depends  on  the  sign  of  the  majority  charge  car¬ 
riers,  the  sign  of  the  electro-optic  coefficient,  and  the  crystal 
orientation. 

This  phenomenon  has  been  reported  in  several  photore¬ 
fractive  media  with  the  most  pronounced  effect  occurring  in 
iron  doped  lithium  niobate  (Fe :  LiNb03)  for  beam  propaga¬ 
tion  and  the  grating  K  vector  along  the  crystal  c  axis  in  a 
focal  plane  geometry.1-4  Although  the  electro-optic  coeffi¬ 
cients  in  Fe:KNb03  are  generally  larger  than  those  of 
Fe:LiNb03,  nonreciprocal  light  transmission  along  the  c  axis 
in  Fe:KNb03  is  less  pronounced  than  with  Fe:LiNb03.  This 
is  because  Fe :  KNb03  does  not  have  a  typically  large  photo¬ 
voltaic  field  which  dominates  beam  coupling  in  most 
Fe:LiNb03  crystals.5  However,  unlike  Fe:LiNb03,  for  this 
geometry  the  effective  electro-optic  coefficient  for 
Fe :  KNb03  increases  dramatically  away  from  the  c  axis  for 
light  polarization  and  grating  K  vector  in  the  crystal  b-c 
plane,  reported  to  reach  a  maximum  at  52.5°  from  the  c 
axis.6 

Experimentally  investigating  the  crystal  angular  depen¬ 
dence  of  self-pumped  nonreciprocal  light  transmission  in 
Fe:KNb03  is  not  trivial.  Photorefractive  Fe:KNb03  is  diffi¬ 
cult  to  grow  reproducibly  and  successive  crystals  grown  us¬ 
ing  the  same  apparatus  can  have  a  wide  variation  in  beam 
coupling  efficiency.  It  is  also  common  for  individual  crystals 
to  have  a  spatially  inhomogeneous  photorefractive  response.7 
The  material  is  fragile  and  sensitive  to  mechanical  and  ther¬ 
mal  stress,  making  it  difficult  to  retain  poling  integrity  during 
cutting  and  polishing  operations.  Re-poling  crystals  cut  sub¬ 
stantially  away  from  the  crystal  c  axis  is  challenging.  Our 
experience  suggests  that  it  is  difficult  to  maintain  experimen¬ 


tal  continuity  across  multiple  samples  in  order  to  make  quan¬ 
titative  measurements.  Ideally  all  comparative  angular  mea¬ 
surements  should  be  made  using  a  single  crystal  to  avoid  the 
large  discontinuities  in  experimental  data  reported  previously 
using  different  crystals.6  This  paper  describes  the  experimen¬ 
tal  measurement  of  self-pumped  nonreciprocal  light  trans¬ 
mission  through  a  uniform  single  crystal  of  Fe :  KNb03  as  a 
function  of  the  crystal  angle  with  respect  to  the  grating  K 
vector  using  a  combination  of  oil  submersion  and  relay  im¬ 
aging.  Using  this  technique  we  show  that  the  optimum  crys¬ 
tal  angle  is  substantially  smaller  than  previously  reported. 
We  also  find  the  beam  coupling  with  crystal  angle  deviates 
significantly  from  published  theory,  suggesting  the  effective 
charge  trap  density  in  Fe:KNb03  is  highly  anisotropic.  Ac¬ 
cording  to  our  measurements,  the  effective  trap  density  for 
the  crystal  a-c  plane  varies  by  almost  two  orders  of  magni¬ 
tude  with  crystal  angle. 

THEORY 

We  define  the  change  in  transmission  due  to  the  presence 
of  a  photorefractive  grating  as  being  a  change  in  the  optical 
density,  AOD,  given  by  AOD  =  log10(/0//p/j),  where  70  is  the 
linear  transmitted  intensity  prior  to  the  formation  of  the  pho¬ 
torefractive  grating  and  IPR  is  the  steady-state  transmitted 
intensity  of  the  pump  beam  once  the  photorefractive  grating 
has  formed.  Using  standard  photorefractive  theory  it  is  pos¬ 
sible  to  estimate  the  optical  gain  coefficient  and  from  this 
calculate  the  expected  change  in  AOD.  Within  the  photore¬ 
fractive  crystal,  standing  intensity  fringes  photoexcite 
charges  which,  after  being  subject  to  thermal  diffusion,  result 
in  a  steady-state  space  charge  field.  Assuming  a  single  ma¬ 
jority  charge  carrier,  Kukhtarev’s  equations  lead  to  the  fol¬ 
lowing  expression  for  the  space-charge  field  Es:5 

g _ Ep  +  iED  +  Epy _ 

5  "  1  +  Ed/Eq  -  i[E0/EQ  +  (Na/Nd)(Epv/Eq)] 

where  E0  is  the  bulk  electric  field,  ED=2irkBT/eA  is  the 
diffusion  field,  EQ=eNA/ Kes  is  the  saturation  field,  EPV  is 
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the  photovoltaic  field,  NA  is  the  acceptor  number  density,  ND 
is  the  total  Fe  number  density  (Fe2++Fe3+),  e  is  the  elec¬ 
tronic  charge,  K=2tt/ A  is  the  spatial  frequency  in  the  me¬ 
dium,  kB  is  Boltzmann’s  constant,  A  is  the  grating  spacing,  T 
is  the  absolute  temperature,  and  es  is  the  dielectric  constant. 
For  Fe:KNb03  the  photovoltaic  field  is  reported  to  be 
negligible8,9  and  so,  in  the  absence  of  an  externally  applied 
field,  Eq.  (1)  reduces  to 


E,=  - 


iEr 


1  +  EJEt 


Q 


(2) 


The  intrinsic  coupling  coefficient,  T,  between  the  two 
counter-propagating  beams  is  given  by 


2tt  ~  . 

T  =  —n3refflm(Es), 

A 


(3) 


where  n  is  the  refractive  index  for  a  given  crystal  angle  and 
beam  polarization,  reff  is  the  scalar  effective  electro-optic 
coefficient,  and  X  is  the  wavelength  of  light  in  vacuum. 

Piezoelectric  and  photoelastic  contributions  arising  from 
the  space-charge  field  can  significantly  affect  beam  coupling 
in  Fe:KNb03.10  The  influence  of  these  contributions  is  de¬ 
scribed  by  an  effective  second  rank  electro-optic  tensor  rft, 
the  components  of  which  are  given  by10,11 

rf  =  A]  A  +  P'ijk,EWmBm  (4) 

in  which 

^ ik~  Cijkl fifth 


E  i  e  kijft  kl^j  ’  (-0 

where  rft  is  the  third  rank  clamped  electro-optic  tensor,  pftft 
is  the  effective  fourth  rank  elasto-optic  tensor,  CEkl  is  the 
fourth  rank  elastic  stiffness  tensor,  ekij  is  the  third  rank  pi¬ 
ezoelectric  tensor,  and  nj  k  l  is  the  unit  K  vector  for  the  space- 
charge  field.  Einstein  summation  over  common  subscripts 
applies  to  these  expressions.  The  scalar  electro-optic  coeffi¬ 
cient  is  then  given  by10 

reff=nP-re//  •%,  (6) 

where  nP  and  ns  are  the  unit  vectors  for  the  pump  and  signal 
beam  polarization  vectors,  respectively.  The  piezoelectric 
and  photoelastic  effects  also  influence  the  dielectric  constant, 
es=  efte{p  where  60  is  the  permittivity  of  free  space  and  eft 
is  the  effective  relative  dielectric  constant  given  by11 

ef  =  t  (7) 

6o 

where  e?-  is  the  clamped  dielectric  tensor. 

The  saturation  field,  Eq ,  is,  therefore,  also  affected  by  the 
piezoelectric  and  photoelastic  effects,  as  well  as  by  the  ef¬ 
fective  trap  density,  NA.  The  influence  this  has  on  beam  cou¬ 
pling  depends  on  the  relative  magnitude  of  the  saturation 
field  with  respect  to  the  diffusion  field  in  Eq.  (2). 

The  coupled  equations  for  the  counter-propagating  beams 


FIG.  1 .  Method  for  determining  T  and  I  Erasure  along  the  crystal 
c  axis. 


db  .  r  I  pis  2  {z-f)IP 

- =  -  alp  -  1 - -  - -  (8) 

dz  Ip  4"  Is  I  Erasure  Zp  (z  f) 

d± s  hh  _2  (z-f)Is 

dz  aS  h  +  h  +  hrasure  Z*  +  (z-/)2’ 

where  IP  is  the  pump  intensity,  Is  is  the  signal  intensity,  a  is 
the  linear  absorption  coefficient,  Zr=tt(oI/X  is  the  Rayleigh 
range  for  the  focusing  lens,  o)0  is  the  He  spot  radius,  X  is  the 
pump  wavelength,  /  is  the  lens  focal  length,  and  z  is  the 
distance  from  the  beam  waist  along  the  beam  axis.  I  Erasure  is 
the  equivalent  intensity  required  to  erase  the  grating  at  the 
same  rate  as  the  sum  of  contributions  from  the  dark  conduc¬ 
tivity  and  any  experimentally  induced  erasure  processes  such 
as  light  scattering  from  defects,  room  light,  etc. 

Equations  (8)  and  (9)  can  be  solved  numerically  using  a 
shoot  and  match  method  with  an  adaptive  Runge-Kutta  inte¬ 
gration  technique  to  estimate  the  AOD.3-5  At  each  iteration, 
the  signal  intensity  is  given  by  the  Fresnel  reflection  of  the 
pump  intensity  at  the  rear  of  the  crystal.  Although  the  spe¬ 
cific  values  of  NA ,  ND,  and  Erasure  are  not  known  for  our 
crystal,  T  and  Erasure  can  be  determined  by  fitting  to  mea¬ 
surements  of  AOD  along  the  crystal  c  axis  for  a  range  of 
input  intensities.  Es  is  then  determined  by  Eq.  (3)  and  we  can 
use  this  value  to  calculate  the  saturation  field,  Eq ,  and  the 
effective  trap  density,  NA.  For  a  given  pump  intensity,  the 
values  of  NA  and  I  Erasure  together  with  the  calculated  value  of 
reff  allow  T  and  AOD  to  be  determined  as  a  function  of  the 
crystal  orientation. 


EXPERIMENT 

T  and  ^Erasure  were  determined  from  measurements  of  the 
AOD  along  the  crystal  c  axis  using  the  arrangement  shown  in 
Fig.  1.  A  5  mm  X  5  mm  X  5  mm  x-y-z  cut  crystal  of 
Fe:KNb03  was  carefully  selected  for  its  uniform  photore- 
fractive  properties.  Unusually,  for  Fe:KNb03,  the  AOD  for 
this  crystal  was  generally  invariant  across  the  whole  pump 
face  of  the  crystal,  indicating  good  volumetric  uniformity  of 
the  photorefractive  properties.  A  532  nm  continuous  wave 
beam  from  a  Coherent  Verdi®  laser  was  focused  at  approxi¬ 
mately  // 40  using  a  +150  mm  focal  length  plano-convex 
lens  with  the  beam  waist  in  the  middle  the  Fe :  KNb03  crys¬ 
tal.  The  beam  polarization  was  parallel  to  the  crystal  a  axis 
and  beam  propagation  was  along  the  positive  c  axis  direc¬ 
tion.  A  large  area  photodiode  was  used  to  measure  the 
steady-state  transmission  for  a  range  of  input  powers.  The 
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FIG.  2.  (Color  online)  Variation  of  AOD  with  input  power  along 
the  c  axis  with  best  fits  for  conventional  and  piezoelectric/ 
photoelastic  theories. 

Fe:KNb03  crystal  had  an  absorption  coefficient  of 
0.74  cm-1  at  532  nm  along  the  crystal  c  axis,  polarized  par¬ 
allel  to  the  a-c  plane.  The  Fe  doping  concentration  in  the 
melt  was  nominally  1000  ppm.  Figure  2  shows  the  variation 
in  AOD  with  input  power  for  this  crystal. 

The  experimental  arrangement  for  measuring  the  crystal 
AOD  angular  dependence  is  shown  in  Fig.  3.  Three  equal 
focal  length  lenses  formed  a  relay-imaged  optical  path 
through  the  apparatus  with  the  middle  of  the  crystal  placed  at 
the  first  focus.  To  allow  higher  internal  beam  angles  to  be 
accessed,  the  crystal  was  submerged  under  Cargille  type  “A” 
immersion  oil  (ft— 1.515  @  532  nm)  contained  within  a 
20  mm  path  glass  cell.  The  glass  cell  was  angled  away  from 
normal  incidence  by  a  few  degrees  to  prevent  the  rear  face  of 
the  cell  from  back  reflecting  a  portion  of  the  pump  beam  into 
the  focal  region  in  the  crystal.  Within  the  glass  cell,  the  crys¬ 
tal  was  mounted  on  a  low  friction  neodymium  magnet.  A 
complementary  magnet  was  positioned  just  below  the  glass 
cell,  fixed  to  a  calibrated  rotation  stage.  Magnetic  coupling 
allowed  the  crystal  to  be  rotated  to  any  desired  angle  [owing 
to  the  aspect  ratio  of  the  crystal  and  to  the  refractive  index 
mismatch  between  the  immersion  oil  and  the  crystal,  internal 
angles  between  32  and  53°  with  respect  to  the  c  axis  could 


f  ,,  f  2f  f 


FIG.  3.  Optical  layout  (top)  and  magnetically  coupled  crystal 
rotation  stage  (bottom).  L=  +  150mm  focal  length  piano  convex 
lenses,  C  =  magnetically  coupled  oil  immersion  cell  for  the  crystal, 
M  =  20%R  mirror,  PD = photodiode,  f- 150  mm  focal  length. 


FIG.  4.  (Color  online)  Conventional  and  piezoelectric/ 
photoelastic  theories  compared  with  a-c  plane  data. 

not  be  accessed  in  our  experiment].  At  the  second  focus  of 
the  relay  imaging  scheme,  a  20%  reflectivity  mirror  gener¬ 
ated  a  counter-propagating  signal  beam  from  the  transmitted 
pump.  This  reflectivity  was  similar  to  the  —16%  reflection 
which  occurs  at  the  air  interface  of  bare  crystals  of 
Fe:KNb03.  Relay  imaging  reduced  the  alignment  sensitivity 
of  the  rear  mirror  and  ensured  a  good  spatial  overlap  of  the 
pump  and  signal  beams.  This  was  especially  important  for 
internal  angles  ^30°  when  beam  fanning  caused  significant 
lateral  distortion  to  the  transmitted  beam.  The  experiment 
was  mounted  on  a  vibration  isolation  table  and,  as  an  addi¬ 
tional  precaution  against  air  currents  and  acoustic  distur¬ 
bances,  the  whole  experiment  was  surrounded  by  foam  rub¬ 
ber  blocks  to  stabilize  the  optical  interference  fringes  within 
the  crystal. 

RESULTS  AND  DISCUSSION 

The  theoretical  fits  to  the  data  for  the  c-axis  power  depen¬ 
dence  shown  in  Fig.  2  give  similar  gain  coefficients  of  16.6 
and  16.1  cm-1  for  conventional  and  piezoelectric/ 
photoelastic  theories,  respectively,  although  the  piezoelectric 
theory  implies  a  greater  intrinsic  erasure  intensity  {I Erasure)  - 
Through  Eq.  (3)  the  values  for  these  gain  coefficients  imply 
space-charge  fields  with  imaginary  components  of  4.70 
X  103  V  cm-1  and  3.86  X  103  V  cm  1  for  conventional  and 
piezoelectric/photoelastic  theory,  respectively,  leading  to  val¬ 
ues  of  Na  of  2.1  X  1016  and  5.0  X  1016  cm-1  from  conven¬ 
tional  and  piezoelectric/photoelastic  theories,  respectively. 
These  are  comparable  with  typical  published  values  for 
Fe :  KNb03  and  are  assumed  to  be  constants  for  the  material. 

Using  the  above  values  for  NA,  Fig.  4  compares  the  ex¬ 
perimental  angular  dependence  of  the  AOD  with  conven¬ 
tional  and  piezoelectric/photoelastic  theories  for  beam  propa¬ 
gation  and  polarization  in  the  crystal  a-c  plane.  It  is  clear  that 
in  neither  case  do  the  theoretical  curves  properly  follow  the 
experimental  results.  There  does  seem  to  be  a  good  correla¬ 
tion  between  conventional  theory  and  experiment  in  Fig.  4 
for  internal  angles  between  zero  and  23°,  but  no  correlation 
exists  for  larger  angles.  There  is  also  reasonable  agreement 
between  piezoelectric/photoelastic  theory  and  experiment  for 
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FIG.  5.  (Color  online)  Conventional  and  piezoelectric/ 
photoelastic  theories  compared  with  b-c  plane  data. 

angles  between  55  and  90°  from  the  c  axis.  Note  that  the 
gain  direction  reverses  at  about  23°  for  conventional  theory 
and  at  about  55°  for  piezoelectric/photoelastic  theory.  The 
crystal  orientation  must  be  reversed  along  the  c  axis  during 
the  experiment  for  angles  greater  than  these  values.  The  data 
presented  allow  for  the  reversal  of  the  crystal  and  absolute 
values  of  AOD  have  been  plotted  in  Fig.  4.  Confirmation 
measurements  were  also  taken  without  reversing  the  crystal 
orientation,  as  a  check  against  orientational  error. 

Using  the  same  respective  values  for  the  space-charge 
field,  Fig.  5  compares  the  experimental  angular  dependence 
of  the  AOD  with  theory  for  beam  propagation  and  polariza¬ 
tion  in  the  crystal  b-c  plane.  For  this  geometry,  conventional 
theory  fails  to  follow  the  experimental  data.  However, 
piezoelectric/photoelastic  theory  gives  a  much  closer  agree¬ 
ment  with  experiment,  but  the  agreement  is  good  only  at 
smaller  angles  and  becomes  progressively  worse  at  larger 
angles.  Note  that  the  peak  beam  coupling  occurs  at  a  beam 
propagation  angle  of  approximately  30-35°  from  the  c  axis, 
which  is  considerably  smaller  than  previously  published 
values.6 

For  both  the  a-c  and  b-c  crystal  planes  the  piezo/ 
photoelastic  theory  gives  a  generally  better  agreement  with 
experiment  than  conventional  diffusion  theory,  although  the 
AOD  (and  hence  the  beam  coupling  gain)  is  suppressed  com¬ 
pared  with  theory  for  beam  angles  away  from  the  crystal  c 
axis.  This  is  particularly  apparent  in  the  a-c  plane  for  angles 
between  15  and  35°.  A  possible  reason  is  that  the  magnitude 
of  the  space  charge  field  is  not  constant  with  crystal  angle, 
although  this  would  be  surprising  in  steady  state  for  charge 
migration  dominated  by  thermal  diffusion.  This  assumption 
changes  if  charge  migration  is  also  driven  by  nonthermal 
mechanisms  such  as  photovoltaic  fields.  The  symmetrical 
beam  fanning  observed  in  our  crystal  at  larger  angles  does 
imply  the  existence  of  a  significant  photovoltaic  field,  in  con¬ 
trast  to  published  literature.8,9  However,  the  photovoltaic 
contribution  to  the  space-charge  field  is  expected  to  have  a 
cosine  dependence  with  crystal  angle  which  does  not  account 
for  the  observed  deviation  from  theories.  Although  the  pres¬ 
ence  of  a  photovoltaic  field  would  adversely  affect  the  grat¬ 
ing  phase  shift  [given  by  the  argument  of  Eq.  (1)  for  the 


FIG.  6.  (Color  online)  Measured  effective  trap  densities  with 
crystal  angle  for  the  a-c  and  b-c  crystal  planes.  Supplementary  data 
are  from  transmission  grating  measurements. 

space-charge  field],  the  net  gain  is  always  enhanced  by  the 
presence  of  a  photovoltaic  field,  described  by  Eqs.  (1)  and 
(3).  This  enhancement  is  greatest  in  the  direction  of  the  pho¬ 
tovoltaic  field,  with  a  cosine  reduction  at  other  angles.  To  test 
for  the  possible  influence  of  the  photovoltaic  field  on  beam 
coupling,  the  phase  shift  of  the  grating  with  respect  to  the 
intensity  fringes  was  measured  by  mounting  the  rear  20% 
reflectivity  mirror  onto  a  piezoelectrically  driven  translator. 
By  translating  the  rear  mirror  along  the  beam  path  at  a  rate 
significantly  faster  than  the  response  time  of  the  experiment, 
a  transient  phase  shift  could  be  introduced  between  the  es¬ 
tablished  photorefractive  grating  and  the  optical  intensity 
fringes.  However,  the  beam  coupling  was  reduced  for  all 
transient  phase  shifts  (both  positive  and  negative  phase  di¬ 
rections),  indicating  an  ideal  90°  phase  shift  existed  in  steady 
state  between  the  index  grating  and  the  interference  fringes 
at  all  crystal  angles.  The  contribution  to  beam  coupling  from 
the  photovoltaic  effect  is,  therefore,  considered  negligible  in 
this  experiment  and  does  not  account  for  the  deviation  of 
beam  coupling  from  theory  at  larger  angles  from  the  c  axis. 

EFFECTIVE  TRAP  DENSITY  MEASUREMENTS 

In  an  effort  to  isolate  the  cause  for  the  theoretical  devia¬ 
tion,  the  measurement  procedure  used  for  Fig.  2  to  obtain  an 
estimate  of  the  effective  trap  density,  NA,  was  repeated  using 
the  apparatus  shown  in  Fig.  3  at  each  measured  angle  for  the 
a-c  and  b-c  crystal  planes.  Assuming  the  published 
values10,11  for  the  electro-optic,  elasto-optic,  elastic  stiffness 
and  piezoelectric  tensors  are  correct,  Fig.  6  shows  the  de¬ 
rived  apparent  effective  trap  density  as  a  function  of  crystal 
angle  for  both  crystal  planes.  There  is  a  surprising  marked 
variation  of  the  effective  trap  density  with  crystal  angle  in 
both  cases,  but  the  apparent  a-c  plane  effective  trap  density 
variation  is  particularly  pronounced,  changing  by  almost  two 
orders  of  magnitude  with  crystal  angle.  Fitted  empirical 
functions  used  to  describe  the  effective  trap  density  varia¬ 
tions  with  crystal  angle  for  the  a-c  and  b-c  planes,  respec¬ 
tively,  are 
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NA(AC)  =  - 9.09  X  lO12#6  +  3.03  X  lO15^  -  3.85  X  lO'V 

+  2.26  X  lO19#3  -  5.67  X  102V  +  3.21  X  1021<9  +  5 
X  1022,  (10) 

NA(bc)  =  -  6  X  10‘V+  1  X  1O2O0+  4.42  X  1022,  (11) 

where  Na^acbq  are  the  effective  trap  densities  in  units  of 
m-3  for  the  a-c  and  b-c  crystal  planes,  respectively,  and  6  is 
the  beam  propagation  angle  in  degrees  with  respect  to  the 
crystal  c  axis. 

We  have  also  considered  the  existence  of  anisotropic  grat¬ 
ings  arising  from  angular  variations  of  the  absorption  coeffi¬ 
cient  and  photoexcitation  cross  section.13-15  However,  the 
variations  in  absorption  coefficient  appear  to  be  insufficient 
to  account  for  our  observations  [a=0.72  cm-1  (along  the 
<2-axis)  to  0.74  cm-1  (along  the  c  axis)  for  the  a-c  plane, 
polarized  parallel  to  the  a-c  plane].  We  have  additionally 
considered  anisotropic  gratings  arising  from  any  anisotropy 
of  the  charge  diffusion  constants  with  crystal  angle.  This 
would  certainly  affect  the  transient  response  of  the  crystal 
but  should  not  influence  the  steady  state  results  upon  which 
our  calculations  are  based. 

Competition  between  electrons  and  holes  can  lead  to  a 
reduction  in  gain,16-18  and  any  crystal  angular  dependence  of 
the  electron  and  hole  diffusion  coefficients  could  lead  to  an 
angular  variation  of  beam  coupling.  However,  our  crystals  do 
not  display  any  change  in  the  sign  of  the  optical  gain  coef¬ 
ficient  with  grating  spacing,  a  key  indicator  of  the  presence 
of  electron-hole  competition.19  The  photorefractive  gratings 
in  our  experiment  decay  with  a  single  time  constant,  indicat¬ 
ing  a  single  majority  charge  species  is  present.  Furthermore, 
competition  between  electrons  and  holes  should  reduce  the 
net  beam  coupling  in  steady  state,  whereas  the  observed 
strong  beam  coupling  supports  the  presence  of  a  single 
charge  species.  For  these  reasons  any  competition  between 
electrons  and  holes  is  thought  to  be  an  unlikely  reason  for  the 
angular  discrepancies  between  theory  and  experiment.  Al¬ 
though  surprising,  anisotropy  of  the  effective  trap  density  is 
therefore  the  most  likely  cause  of  the  deviation  from  previ¬ 
ously  published  theories. 

We  have  made  a  concerted  effort  to  independently  verify 
the  unprecedented  variation  in  the  effective  trap  density  us¬ 
ing  the  standard  method  of  measuring  the  small  signal  gain 
as  a  function  of  the  grating  spacing  using  a  combination  of 
transmission  and  reflection  grating  measurements.19  The 
method  relies  on  identifying  the  optimum  grating  spacing  to 
determine  the  effective  trap  density  of  the  material.  By  re¬ 
peating  these  measurements  over  a  range  of  crystal  angles, 
the  effective  trap  density  as  a  function  of  the  crystal  angle 
and  crystal  plane  can  be  obtained.  Owing  to  the  tensoral 
nature  of  the  electro-optical,  piezoelectric,  and  photoelastic 
coefficients,  it  is  not  possible  to  obtain  an  analytic  expression 
for  the  optimum  grating  spacing  for  each  crystal  angle.  Our 
approach  has,  therefore,  been  to  adjust  the  scalar  effective 
trap  density  for  each  given  crystal  angle  until  the  theoretical 
optimum  grating  spacing  coincides  with  the  experimental  ob¬ 
servations  for  that  crystal  angle.  By  repeating  this  process  for 
each  crystal  angle  and  crystal  plane,  the  effective  trap  density 


can  be  extracted  as  a  function  of  the  crystal  orientation. 

This  method,  although  simple  in  concept,  is  experimen¬ 
tally  extremely  challenging  in  our  situation.  As  the  crystal  is 
rotated,  the  relative  angles  of  incidence  of  the  pump  and 
signal  beams  become  asymmetrical,  leading  to  rotation  of  the 
grating  K  vector  with  respect  to  the  crystal  axes.  The  me¬ 
chanical  angle  of  the  crystal,  therefore,  required  adjustment 
for  every  change  of  the  beam  intersection  angles  in  order  to 
maintain  a  constant  grating  vector  orientation.  It  was  also 
very  important  to  maintain  a  good  spatial  overlap  between 
the  pump  and  signal  beams  for  all  beam  intersection  angles 
and  crystal  orientations.  This  demanded  making  the  pump 
beam  larger  in  diameter  than  the  signal  beam.  The  relatively 
large  diameter  of  the  pump  beam  made  it  necessary  to  rotate 
the  crystal  towards  the  pump  in  order  to  avoid  clipping  of  the 
pump  beam  by  the  crystal  edges.  The  large  diameter  pump 
beam  and  large  gain  coefficients  promoted  excessive  compe¬ 
tition  between  signal  amplification  and  parasitic  photorefrac¬ 
tive  coupling  arising  from  scattered  pump  light,  leading  to 
strong  beam  fanning  and  spurious  beams.  It  proved  impos¬ 
sible  to  make  meaningful  measurements  at  crystal  angles  be¬ 
yond  25°  in  the  a-c  plane  or  15°  in  the  b-c  plane.  We  at¬ 
tempted  to  reduce  the  overall  gain  using  an  external 
incoherent  light  source  to  compete  with  the  grating  forma¬ 
tion,  but  this  was  ineffective  and  we  could  not  prevent  com¬ 
petition  from  parasitic  photorefractive  processes.  These  are 
the  primary  reasons  we  used  the  method  shown  in  Fig.  2  to 
obtain  trap  density  information,  which  we  believe  to  be  more 
reliable  than  the  results  obtained  from  the  transmission  grat¬ 
ing  method  described  above.  However,  the  transmission  grat¬ 
ing  trap  density  measurements  did  show  qualitative  agree¬ 
ment  with  our  reflection  grating  observations.  The 
transmission  grating  results  are  shown  as  supplementary  data 
in  Fig.  6.  There  is  clear  evidence  that  the  b-c  plane  effective 
trap  density  remains  approximately  constant  with  crystal  ori¬ 
entation,  whereas  the  a-c  plane  effective  trap  density  appar¬ 
ently  falls  sharply  as  the  grating  K  vector  is  rotated  away 
from  the  crystal  c  axis,  supporting  our  earlier  observations 
from  self-pumped  reflection  gratings. 

The  effective  trap  density  data  from  Fig.  6  (reflection 
grating  values)  can  be  used  to  produce  a  revised  theoretical 
estimate  of  the  AOD  with  crystal  angle  by  using  the  fitted 
empirical  functions  for  the  effective  trap  density  variations 
for  the  a-c  plane  [Eq.  (10)]  and  for  the  b-c  plane  [Eq.  (11)]  in 
place  of  the  previously  fixed  value  of  NA.  Using  this  ap¬ 
proach,  Fig.  7  shows  the  calculated  AOD  for  both  crystal 
planes  and  gives  strong  agreement  with  experiment.  Note 
that  the  trap  density  data  was  obtained  from  a  different  ex¬ 
periment  and  data  set  than  the  original  AOD  measurements, 
so  that  the  good  agreement  between  theory  and  experiment  is 
not  due  to  a  circular  argument  from  using  common  data. 

The  apparent  presence  of  an  anisotropic  trap  density  is 
rather  surprising  and  has  important  implications  for  semicon¬ 
ductors  in  general,  not  just  photorefractive  materials.  An  im¬ 
plausible  explanation  is  that  the  trap  sites  are  somehow  sen¬ 
sitive  to  the  polarization  state  of  the  original  photon  used  to 
excite  the  electron  (or  hole)  into  the  conduction  (or  valence) 
band.  Each  charge  would  need  to  retain  information  about 
the  original  exciting  photon.  It  is  widely  assumed  that  this 
information  is  lost  once  an  electron  becomes  delocalized  in 
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FIG.  7.  (Color  online)  AOD  variation  with  crystal  angle  using 
the  angular  effective  trap  density  measurements  for  the  a-c  and  b-c 
crystal  planes. 

the  conduction  band.  If  information  of  the  original  exciting 
photon  is  not  lost,  as  has  been  assumed,  then  this  could  ac¬ 
count  for  the  observed  trap  density  anisotropy.  However,  we 
believe  this  to  be  unlikely  because  the  electron  spin  relax¬ 
ation  times  are  orders  of  magnitude  faster  than  the  electron 
trapping  times. 

As  an  alternative  explanation  for  the  apparent  anisotropy, 
we  propose  that  the  effective  trap  sites  within  KNb03  may 
occupy  a  complex  “layered”  structure.  Growth  of  KNb03 
does  not  occur  smoothly  with  a  constant  growth  rate  and 
instead  proceeds  in  a  stochastic  manner  adding  a  plethora  of 
microscopically  thin  layers  of  new  material  to  the  existing 


bulk  single  crystal.  It  has  been  observed  that  impurities  tend 
to  segregate  at  the  boundaries  between  these  structures.20  It 
is  reasonable  to  assume  that  the  added  impurities  (Fe)  form  a 
stratified  concentration  of  traps  (Fe3+  and  intrinsic  defects)  at 
these  boundaries.  If  the  mean  distance  between  adjacent  lay¬ 
ers  varies  with  crystal  angle,  the  effective  trap  density  could 
become  highly  anisotropic  when  the  mean  layer  separation 
becomes  comparable  to,  or  greater  than,  the  optical  fringe 
separation.  For  gratings  with  K  vectors  predominantly  paral¬ 
lel  to  the  layers,  the  average  trap  density  would  appear  ap¬ 
proximately  constant.  Whereas  for  gratings  with  K  vectors 
generally  orthogonal  to  the  layers,  some  fringes  would  expe¬ 
rience  very  low  effective  trap  densities  while  others  would 
have  higher  trap  densities.  For  this  situation,  the  local  space- 
charge  field  may  saturate  in  the  low  trap  density  regions, 
leading  to  an  overall  reduction  in  beam  coupling.  We  note, 
however,  that  the  proposed  layering  of  the  trap  structure 
must  be  irregular  in  nature;  a  regular  spatial  frequency  would 
lead  to  an  intrinsic  trap  grating  in  the  crystal,  leading  to 
competition  with  the  photorefractive  grating  when  the  re¬ 
spective  K  vectors  become  parallel.  An  intrinsic  trap  grating 
should  be  observable  using  incoherent  light,  but  we  did  not 
find  any  such  gratings  in  our  crystal. 

CONCLUSIONS 

The  dependence  of  self-pumped  photorefractive  two- 
beam  coupling  in  Fe:KNb03  with  crystal  angle  has  been 
characterized  for  the  a-c  and  b-c  crystal  planes  at  532  nm. 
Existing  piezoelectric/photoelastic  theory  gives  better  agree¬ 
ment  with  experiment  than  conventional  diffusion  theory,  but 
significant  deviations  occur  at  large  angles  from  the  c  axis, 
especially  for  the  a-c  crystal  plane.  We  find  evidence  to  sug¬ 
gest  that  the  effective  charge  trap  density  in  Fe :  KNb03  may 
be  highly  anisotropic  with  almost  two  orders  of  magnitude 
variation  with  crystal  angle  in  the  a-c  plane,  leading  to  a 
tensoral,  rather  than  scalar,  description  of  the  effective  trap 
density  in  Fe:KNb03.  We  propose  that  the  apparent  trap 
density  anisotropy  could  be  due  to  a  concentration  of  traps  at 
the  interface  between  adjacent  irregular  growth  layers  during 
crystal  production. 
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